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Abstract - 

Langmuir monolayers are modeled as systems of short chains, which are 
confined to a planar surface at one end, but free to move within the plane. The 
phase behavior is calculated in a mean field approximation, which combines 
the self consistent field method with elements of classical density functional 
theory. It is shown that phases with tilt order are unstable in systems of 
stiff chains, but can be stabilized by chain conformational entropy in systems 
of sufficiently flexible chains. The chain entropy is also responsible for the 
appearance of an additional untilted phase, the liquid expanded phase. The 
region of stability of the different phases is discussed, and their microscopic 
structure is analyzed in some detail. 
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1 Introduction. 



Monolayers of amphiphilic molecules have been studied for many years for 
practical and fundamental reasons. Placed on a solid substrate, they build 
Langmuir-Blodgett films, which have important technical applications, e.g., 
in thin film technology JT|. Monolayers of lipids on water are of biological 
interest, since lipid bilayers - consisting of two weakly coupled monolayers - 
are essential ingredients of biological membranes M. 

The phase diagram of Langmuir monolayers (monolayers adsorbed at the air 
water interface) at low surface coverage is qualitatively similar for long chain 
fatty acids, alcohols and lipids (Figure 1) |5|, ||. It's most remarkable feature 
is the presence of two distinct fluid-fluid coexistence regions at intermediate 
temperatures: the familiar "gas - liquid" transition at low surface densities, 
and an additional transition from a "liquid expanded" (LE) phase to a "liquid 
condensed" state at higher densities. The latter is indeed a transition between 
two fluid states, as evidenced by the experimental observation that positional 
correlations in the condensed phase decay exponentially ||. It is the mono- 
layer equivalent of the "main" transition in bilayers, which is interesting from 
a biological point of view, because it is found at temperatures often close to 
the body temperature (41.5 °C in DPPC) At even higher surface coverage, 
monolayers can display a rich spectrum of condensed phases, which differ from 
each other in positional order, tilt order, and orientational order of the back- 
bones of the chains 0. In this work, we shall discuss the condensed phases 
which can coexist with the expanded phase, i.e. the high temperature untilted 
phase (LS) and the low temperature tilted phase (L 2 , see Figure 1). 
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The nature of the transition between liquid expanded and condensed phases 
has been discussed over many years. In an earlier paper ||, we have presented 
self consistent field calculations of a "minimal" model for Langmuir monolay- 
ers, where the amphiphilic molecules were modeled as semiflexible chains with 
one end grafted to a planar surface. We have shown that two ingredients are 
needed to bring about coexistence between two liquid states: The chain flex- 
ibility, which stabilizes the expanded phase, and the chain anisotropy, which 
dominates the liquid condensed state. The transition is driven by the interplay 
between the entropy of chain disorder and the energy associated with collective 
chain alignment. The latter may result from simple packing effects, or from 
additional (e.g., dipolar) anisotropic interactions between chain segments. 

The model hence successfully reproduced the LE and the LS phase, yet 
it seemed to fail to display stable phases with collectively tilted chains. In- 
dications for tilt order were only seen in the unstable regions of two phase 
coexistence. In that respect, the observed phase behavior was similar to that 
of grafted rigid rod systems. Grafted rods with fixed grafting points may show 



tilt order in a region of surface coverage |10|, [TT . However, the surface energy 
per chain is higher in the tilted region than in the untilted region. When 
the rods are given translational degrees of freedom, the tilting transition is 



therefore replaced by phase separation ]I2, 13 ]. 

According to a common picture, tilt order in Langmuir monolayers results 
from a mismatch between head group and tail segment size. The larger area of 
the head group constrains the coverage of the condensed phase and stabilizes 
surface coverage regions with tilt order. This mechanism is doubtless the 
driving force for tilt order in many cases, but it is certainly not the only one. 
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For example, it hardly explains the experimental observation of tilt order in 
monolayers of triple chain phospholipids ||14|| . Another potential cause for 
tilt is related to the internal structure of the chains: When the chains are 
tilted, monomers can "hook" into each other, and thus pack more effectively. 
Presumably, this is responsible for the presence of tilt order in Monte Carlo 
simulations of endgrafted bead-spring-chains [|15|, 16 . Tilt order may also be 
induced by attractive interactions between the chains and the bare surface 
IU| , O . For hydrophobic ((CH 2 ) n ) chains on a water surface, that seems 



however less likely. 

All these tilting mechanisms do not operate in the minimal model of refer- 
ence 0. Therefore one would not expect to find tilt order there, unless the 
model is extended in a suitable way. Yet we shall show that the conformational 
degrees of freedom of the chains generate a new mechanism for the stabilization 
of tilted states: As we have discussed above, giving the chains some flexibility 
brings a new phase into existence, the liquid expanded phase. On making the 
chains more and more flexible, the condensed phase is affected too: The gain of 
conformational entropy at lower surface densities compensates in part the loss 
of surface energy. As a result, the region of stability of the condensed phase 
is extended. Provided the chains are sufficiently anisotropic, the coverage at 
coexistence becomes low enough to support collective tilt. 

The influence of conformational chain disorder on tilt in fatty acid mono- 



layers has received some interest recently |18|, [19| . It has been argued that in 
tilted phases, an increase in the number of gauche defects in the chains reduces 
the tilt angle at the same area coverage. The present work discusses an an- 
tipodal, although related effect: Chain disorder stabilizes homogeneous tilted 
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phases at molecular areas, where ordered, straight, chains phase separate into 
two untilted phases. 

The purpose of this work is two fold: To explore the possibilities for tilt order 
within the minimal model, and to establish a complete phase diagram in terms 
of the variables stiffness and chain anisotropy. The parameter region, in which 
tilted phases are stable, will be determined, as well as the parameter region, 
in which a liquid expanded and a liquid condensed phase can coexist. Where 
those two regions overlap, one finds a phase diagram which is very similar to 
the one sketched in Figure 1. The paper is organized as follows. The model 
and the self consistent field method are described in the next section. A variant 
of the model || is used, which allows among other for a more detailed study 
of chain defects. Section three presents the predictions of the self consistent 
field theory first for flexible chains, then for stiff chains. The properties of the 
different phases are discussed in some detail (density profiles, nematic order, 
chain defects), and an overview over the phase behavior is given. The results 
are summarized in section four. 

2 The Model. 

A schematic picture of the model is shown in Figure 2. The amphiphilic 
molecules are modeled as chains containing n rod-like tail segments of length 
Iq and diameter A , and one head segment, which is confined to a planar surface 
at z = 0. They are subject to three different types of potentials: 

• External potentials, which confine the head segment at z < and the tail 

segments at z > 

5 



• A bending potential, which favors parallel alignment of adjacent segments. 

• The interactions between segments. Tail segments are anisotropic, have 
a repulsive hard core, and attract each other at larger distances. Head 
segment interactions are isotropic and purely repulsive. 

The external potentials h e ^ 1 (head segments) and h\ xt (tail segments) are 
taken to be simply harmonic. 

K\r) _ j z<0 K\r) _ j k t z 2 z<0 

k B T ~ \ hz 2 z>0 and k B T ~ \ z > ' ^ 

where k B is the Boltzmann constant and T the temperature. 

The choice of the bending potential is guided by the idea that tail segments 

in the model chain correspond to two CH 2 groups each in a hydrocarbon chain. 

A molecule in an all trans conformation is then represented by a completely 

stretched model chain with bending angles theta = 0. A conformation with 

one gauche kink is represented by a model chain, which has one bending angle 

6 = (7r/3) or cos8 = (1/2). In view of these considerations, the bending 

potential is given the the form U(Q) j \k B T) = uU(9), where u is an adjustable 

stiffness parameter, and 

U{6) = 25x + 34a; 2 - 400x 3 + 480x 4 with x = 1 - cos 9. (2) 

The function U(8) is plotted in Figure 2. It has a minimum at cos(#) = 1/2, 
and takes the value Uq = 1 there. The relative potential barrier U m /Uo ~ 4 
has approximately the same height as the energy barrier from trans to gauche 
in popular polyethylene models (e.g. by Rigby and Roe 01). Moreover, the 
thermal average of cos(#) in a free model chain at u = 1, (cos(#)) ~ 0.7, is in 
rough agreement with the value obtained for polyethylene at k B T = E g , where 



Eg is the energy of a gauche defect (calculated in the RIS scheme, [^1 



The mapping of carbon groups on chain segments should not be taken 
too literally, since the model is so simple compared to a real hydrocarbon 
chain. However, the choice of a bending potential with two minima such as @ 
has the advantage, that it allows to define chain defects, and to study defect 
distributions. Note that the energy of gauche defects is of order 300-fT, i.e., 
room temperature, in units of the Boltzmann constant ks- Hence an analysis 
of their distribution can be instructive, especially in short chains. For most 
other purposes, a simple harmonic potential such as has been used in Ref. M 
is entirely sufficient, and yields qualitatively the same results. 

The interaction between segments are introduced in terms of a functional 
^[{Phifi w),pt(r, w)}] of the center of mass densities of head (ph) and tail (p t ) 
segments with orientation w (\w\ = 1) at position f. The functional includes 
short range repulsive hard core potentials as well as longer range attractive 
interaction tails. 

In an exact treatment of the above model, one has to perform ensemble av- 
erages over all possible configurations of chains, and the corresponding center 
of mass densities. In this work, we will resort to a local mean field approxima- 
tion. The densities ph,t(r, w) are replaced by their ensemble averages, and T 
is taken to be a functional of average densities. Single segments interact with 
others via average fields 



The effect of local density fluctuations is neglected. 

This approximation has a number of important implications. First, corre- 
lations between different chains are neglected. The problem therefore reduces 
to calculating the partition function and the density distribution of a sin- 




ma! 



(3) 
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gle noninteracting chain (random walk) in the inhomogeneous external fields 
hh }t (r,w) = h™£ + h'jfj:, which have to be determined self consistently using 
eqn @ (cf. [p2 , |23||). Note that correlations within a chain are still present 



due to the chain connectivity. Second, the non integrable hard core interac- 
tions require special treatment. We will choose a common approach in density 
functional theories |24, 25 1, which is to expand around a reference system of 



purely repulsive segments. Third, the mean field approximation does not cap- 
ture the fact that stiff chains are always anisotropic, even if the constituting 
segments are not. Effective anisotropic interactions result, e.g., from packing 
effects. Within the mean field approach, they have to be introduced explicitly 
in terms of an effective segment anisotropy. 

Since the segments are extended objects, the center of mass density p cor- 
responds to a segment mass density 

Ph,t(r,w) = dr'K h j(r-r\w)p h , t (r',w). (4) 



The function K(r, w) reflects the shape of a segment with orientation w. Since 
the segments are fairly compact, the orientation dependence of the shape func- 
tion K(r, w) can be neglected, and it is reasonably well approximated by a 
simple step function. 

m ^ = | 1/Mo) W<V2, (* 2 + y 2 ) < A /n (5) 

We shall also need the total density 

P(f) = J- / d ™ \Pt(r, w) + p h (r, w)}, (6) 

where the integral / dw is performed over the full solid angle Air. With these 
definitions, we are able to formulate a concrete Ansatz for the density func- 



tional F. We use a local density approximation, i.e., the functional F is given 
as the integral over a free energy density function. 

k^f' F = I rff {^°[ p( ^ + 32^2 / / d ™ d ™'Pt{r,w)p t {r,w')[V(w-w')-e)} (7) 
The first term describes a reference system of identical segments with isotropic 
hard core interactions. The free energy density fo[p] is derived from the hy- 
pothetical equation of state of a dense melt of such "ideal" chain segments: 
Being part of a chain, the segments have no translational degrees of freedom, 
their equation of state has no ideal gas contribution. Furthermore, segments 
are connected to others at both ends, therefore they mainly interact within a 
plane perpendicular to themselves. Hence we assume that their equation of 
state is reasonably well approximated by the equation of state for hard disks 



PS|, from which the ideal gas term has been subtracted. 



n(p) = p( 7I ^ - 1) (8) 

with the reduced pressure II = p/k B T, the density p and the packing frac- 
tion 7] = p ■ AqIq. From this one can derive the free energy density using 
d(fo/ P )/dp = n/ P 2 . 

/o[p]=p{ T ^--log(l-7?)} (9) 

1 — 7] 

The second term in eqn (^) accounts for the anisotropic and attractive 
interactions between tail segments perturbatively, up to the leading order in 
the densities. The attractive part of the interaction is absorbed in a single 
parameter e. The anisotropic part of the interaction is described by an even 
function V(x) = V(—x) and can be expanded in Legendre polynomials. 

00 2/ 4- 1 

V(x) = V(-x) = £ —^P^vl (10) 

l=2.4.- 4?F 



We shall neglect all contributions except for the lowest, v = — v-i- It should 
be emphasized again that the anisotropy parameter, v, cannot necessarily be 
traced back to an actual anisotropy of single free segments. It is an effective 
parameter, which has to be introduced in a mean field theory in order to 
include effects of the chain anisotropy. Thus it has to be identified with an 
effective anisotropy per segment, rather than with the anisotropy of a, segment. 

We complete the definition of the model by specifying the parameters 
k h = k t = 20/3 / 2 , e = 40 Zg, A = 2.01 (see ||). This choice is motivated 
as follows: The parameters kh and k t can be chosen arbitrarily, provided 
they are large enough to ensure the confinement of the heads at the sur- 
face, and of the chains above the surface. The strength of the attractive 
interaction, e, determines the density within a hydrophobic layer, and af- 
fects the jump in the surface coverage at first order transitions. As shown 
in reference |J, e has not much qualitative influence on the phase behav- 
ior, therefore it is not varied systematically here. In a virial expansion, e is 
given by the integral over the Mayer-f-function of the attractive interaction, 
e = J(exp[—v a ttr.( r *)/^BT] — l)dr. The parameters e and the effective chain 
diameter A were chosen such that they are compatible with the size and po- 
tentials of alkane chains, if one maps two (CH2) n groups on one model segment, 
with alkane potentials taken from ref. f27| . 

All calculations were done with chains of tail length n = 7. Free model 
parameters, which were systematically varied, are the stiffness parameter u 
and the anisotropy parameter v , hereafter given in units of Zq. We shall com- 
ment briefly on their connection with interaction parameters in other systems, 
e.g., simulation models. In a simulation, the effective stiffness u can be es- 
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timated from matching the thermal average (cos 8) for the angle 9 between 
adjacent bonds, in a dense melt of free chains, with the average obtained for a 
random walk of rods with the bending potential uU (fj). The least accessible 
parameter is the effective anisotropy parameter v. As noted earlier, it's origin 
is mostly due to packing effects. A lower bound can be calculated from the 
excluded covolume (the second virial coefficient) of two stretched chains of the 
persistence length, divided by the number of segments. In the present model, 
at u ~ 1—2, one gets v ~ IOZq. Such a calculation however neglects the 
anisotropy in the attractive interaction. Moreover, the segment density in the 
hydrophobic layer is very high (see Figure 10), such that higher order virial 
coefficients come heavily into play. Hence the resulting effective anisotropy 
will be much higher. In a simulation, v can be determined from the analysis 
of orientation correlations between segments in a melt of free chains. 

The procedure used to solve the problem is similar to the Scheutj ens- Fleer 



method for lattice models of polymers at surfaces p2[ . One defines recursively 
the end segment distributions (i < n) 

Wi(r,w) = I ' dtv'W i . l {f\wy^ h ^- u ^^/ kBT 
4tt J 

f' — f — [w + w'] (11) 

'\J-ht(f',w)-U(w-w'))/k B T 



Wi(r,w) = i- J 'dw'W i+1 {f',w')e 



f' = f + -[w + w') (12) 

with Wq{t,w) = exp(—h h (r)/kBT) and W n (r, w) = 1. We consider a homo- 
geneous monolayer of ./V chains, which occupy each an area per molecule A. 
Hence we have translational invariance on the xy plane, and the single chain 
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partition function is given by 

1 



dzdwWi(z,w)Wi(z,w), (13) 



u 4tt/ 

which is independent of i. The center of mass density of the zth segment can 
be calculated via 

-/ 1 Wj{z,w)Wi(z,w) 

p l (z,w) = — (14) 

and the total free energy per chain (with the de Broglie wavelength X B ) 

= _ teg Zo _ l g(A/A|) - 1 (15) 
f dfo 1 /" 

J dz {(f [p\ ~ P-^) ~ 7^~i J dwdw'p t (z,w)p t (z,w')(V(ww') -e)}. 

The chemical potential u, i.e., the free energy gain on adding one chain is 
therefore given by 

1 dF n Z A 

log— J-. (16) 



k B T dN k B T to A| 

In the grand canonical ensemble, this leads to the Gibbs free energy per surface 

area g 

g 1. F n 



k B T A y Nk B T k B T^' ^ 
Unless stated otherwise, the free energy and the chemical potential will be 

given in units of k B T and shifted by log(A^) — 1 in the following. 

In practice, it is useful to expand functions of orientation w in spherical 
harmonics. Moments up to I = 10 were taken into account, i.e., 121 functions, 
a number which proved sufficient. The z-direction was discretized in steps of 
Zo/5. The mean field equations were solved iteratively, using the Legendre co- 
efficients of the fields h™f (r, w) as iteration variables. The iteration procedure 
combines a method proposed by Ng and simple mixing: Let the vector x n 



be the nth guess of the set of iteration variables, f n the fields calculated from 
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there, and d n = f n — x n the remaining deviation. Following Ng, we define the 
matrix Uij = (d n — d n -i) ■ (d n — d n ^j) and the vector Vj = (d n — d n ^j) • d n , 
where the dimension of U and V, i max = j ma x, is arbitrary (2 to 5 in this 
work). We then invert U, determine the coefficients At = Vj, and calcu- 
late x£ = x n + Y,iA(x n ^i-x n ) and f A = f n +T,iA(fn-i-fn)- In the iteration 
procedure suggested by Ng, the (n + l)th guess of x is given by x n+ i = f A . 
Unfortunately, this method does not converge for the present problem. Good 
results were however obtained with the prescription x n+ \ = x A + A(/^ — x^), 
with A ranging between 0.1 and 0.2. A relative accuracy of 10 -8 was usually 
reached within less than 100 iteration steps. The iteratively obtained solutions 
for fixed surface coverage were usually unique, unless metastable states (e.g., 
tilted states) existed. In that case, the solution with the lowest free energy 
QT5D was selected. 

3 Results. 
3.1 Stiff Chains. 

Figure 3 shows a free energy curve in a system of relatively stiff chains (u = 2, 
v = 13.7). On increasing the molecular area, the free energy exhibits two min- 
ima and then rises. As the area tends to infinity, not shown, it diverges nega- 
tively, following the ideal gas term — \og(A/X^). Hence the Maxwell enveloping 
function has a negative slope, which guarantees the mechanical stability of the 
system: The spreading pressure U/ksT = —dF/dAN" 1 is always positive. 
The pressure in the gas phase is however very low, in the Maxwell construc- 
tion the common tangent with a coexisting gas phase is practically horizontal. 
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Figure 3 illustrates the situation where one has two distinct regions of phase 
separation, first between a condensed phase (LS) and an expanded phase (LE), 
and then between an expanded phase and the gas phase (G). The fact that 
there is phase separation can be inferred from the Maxwell construction, which 
does not follow the free energy curve in those two regimes, and from the obser- 
vation that the Gibbs free energy (|i~7|) is not a unique function of the chemical 
potential (Figure 3, inset). 

If one decreases the chain stiffness u or the chain anisotropy v, the coexisting 
expanded and condensed phases merge into one at a critical point (Figure 4a 
and b). This point is difficult to locate from just looking at the free energy 
curves, but can be identified via the inspection of the Gibbs free energy as a 
function of the chemical potential. On increasing u or v, on the other hand, 
the free energy minimum belonging to the condensed phase decreases relative 
to the other minimum. A triple point is encountered, beyond which the liquid 
expanded state is metastable, and the condensed phase coexists with the gas 
phase. A state with collective tilt emerges in the unstable surface coverage 
region. Tilt order thus occurs in a system of fixed grafted chains, but is 
replaced by phase separation when the chains are allowed to move. 

So far, our results essentially confirm and complete the results reported in 
reference ||. Systems of relatively stiff chains qualitatively show the same be- 
havior as found earlier in a somewhat different model. Hence we shall not dis- 
cuss this regime in more detail. The phase diagram in the plane of anisotropy 
vs. molecular area at chain stiffness u = 2 is shown in Figure 5. 
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3.2 Flexible Chains. 

The transition between the condensed phase and the expanded phase is gov- 
erned by the interplay of chain flexibility and chain anisotropy. In systems of 
more flexible chains, one recovers two phase coexistence if the higher confor- 
mational entropy is compensated by higher effective segment anisotropy. 

Free energy curves for the set of parameters u = 1, v = 20.3, are shown in 
figure 6. As in figure 3, there are two successive first order transitions between 
fluid phases, passing from the gas phase (G) via a liquid expanded phase 
(LE) to a liquid condensed phase (L2). Contrary to the case of stiffer chains, 
however, the coexisting condensed phase is tilted. Upon further compression 
of the monolayer, an additional continuous transition to an untilted state (LS) 
takes place. 

The tilt order can be measured in terms of the in-plane alignment of seg- 
ments, d|| = ,J {vj x ) 2 + (w y ) 2 . The fact that d\\ 7^ implies that the symmetry 
in the xy plane is broken. Figure 7 demonstrates that the "tilted state" indeed 
displays this kind of azimuthal order. In systems of fixed grafted chains, the 
tilted state is stable in a coverage interval, bounded by a continuous transition 
at high coverage (marked I in figs 6-8) and by a first order transition to the 
untilted state at low coverage (marked II). When the chains are given lateral 
mobility, this second transition disappears in the coexistence region of the L2 
and the LE phase. 

Further insight can be gained from the inspection of the nematic order in 
the system. The relevant quantity here is the traceless ordering matrix |H| 
S = (3wiWj — Sij)/2. It has the eigenvalues {S, — (S — r])/2,—(S + ??)/2}, 
with the nematic order parameter S and the biaxiality rj. The nematic order 
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S is always nonzero, since the chains are always aligned to some extent in 
the direction perpendicular to the surface. As the molecular area increases, it 
decreases monotonically in both the tilted and the untilted phase, yet it stays 
higher in the tilted phase. At the first order transition (II), S jumps from 0.37 
in the tilted phase to 0.21 in the untilted phase. The value of S in the tilted 
state is thus comparable to its value in the nematic phase of liquid crystals, 
right at the transition to the isotropic phase (S = 0.43 in the Maier-Saupe 
model |2IJ). The in-plane symmetry breaking is reflected by the behavior of 
the biaxiality, which is nonzero only in the tilted state. 

From these results the nature of the tilting transitions in the system can 
be inferred. The discontinuous low coverage transition (II) is associated with 
ordering/disordering of single segments. It is thus essentially a a nematic- 
isotropic transition, analogous to those found in liquid crystals. The contin- 
uous high coverage transition (I), on the other hand, results from in-plane 
ordering/disordering of whole chains. The surface induces an orientation di- 
rection, hence the transition is of XY type |Kj . The two types of transitions 
are illustrated in figure 9. 

The structure of the monolayer shall be analyzed in some more detail. 
The density profiles in the three phases do not differ remarkably from each 
other. Examples are shown in figure 10. The total segment density is constant 
throughout the layer and independent of the surface area per chain or the tilt 
order. It is also independent of the chain stiffness and chain anisotropy, and 
only determined by the interaction parameter e (not shown). Compression 
results in thickening of the monolayer. 

The distribution of bending angles 8, shown in figure 11, is more interesting. 
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It has by construction two maxima, one at the bending angle 8 = and one 
at cos# = 0.5. The area under the second maximum gives the concentration 
of conformational (gauche) defects in the chains. As demonstrated in the in- 
set, the distribution for the outermost angle, the angle between the last two 
segments, is always the same up to the molecular areas which were consid- 
ered. The main graph shows the deviations from this distribution for the inner 
angles. In the expanded phase, chains have more defects in the middle than 
at the ends, i.e. they are more disordered there. In the condensed phases, 
in contrast, the conformational order is highest in the middle. This result is 



in agreement with molecular dynamics simulations [31, 32] and other model 



calculations 33 



We close this section with the discussion of the phase diagram at chain 
stiffness u = 1 (figure 12). At low chain anisotropy v, there is only one single 
untilted liquid phase, which coexists with the gas phase. A tilted phase L 2 
emerges at a tricritical point, v = 20.1, and separates two untilted liquid 
regions, the expanded (LE) and the condensed phase (LS). The transition 
between the L 2 and LS state is continuous at lower values of the anisotropy v, 
and replaced by phase separation at the tricritical point v = 22.7. We note that 
the tilting transition in monolayers of chains with fixed homogeneous grafting 
density remains continuous. The tilt order parameter vanishes continuously at 
a critical line, which is however hidden in the coexistence region if the chains 
are mobile. Beyond the triple point, where the liquid expanded phase becomes 
metastable (v = 20.5), the region of stability of the tilted phase narrows down 
and finally disappears at v = 34. 
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3.3 Phase behavior. 

We have seen that systems of stiff chains exhibit fluid fluid coexistence of 
two untilted liquid phases, whereas in systems of flexible chains, the liquid 
expanded phase coexists with a tilted condensed phase. In an intermediate 
range of stiffness, one can find both. An example is the phase diagram for 
chains of stiffness u = 1.5, shown in figure 13. Phase separation between an 
expanded phase and an untilted condensed phase sets in at the critical point 
v = 16.8. A tilted phase emerges at v — 17.1 in the coexistence region between 
the expanded and the untilted phase. The liquid expanded phase ceases to be 
stable at the triple point v = 17 A. The L 2 phase and the LS phase are 
separated by a narrow coexistence region; in systems of slightly more flexible 
chains, u < 1.45, the transition can also be continuous in a window of v (see 
figure 14). 

Figure 14 summarizes the phase behavior for chain stiffnesses ranging be- 
tween between u = 1 and u = 2. It shows a projection of the three dimensional 
phase diagram in the (A, u, v) volume into the (w, v) plane. The shaded area 
designates the region where a tilted L 2 state is stable. The transition from 
this phase to the untilted LS phase is continuous in the light shaded area, 
and first order in the dark shaded region (z.e., the two phases phase separate). 
The light shaded area is thus bounded by two lines of tricritical points. The 
liquid expanded phase (LE) is stable in the hatched area, which is bounded 
by a tricritical or critical line, and a triple line. At very large chain stiffness, 
u ~ 3.5, these two lines merge and disappear (not shown). Hence we recover 
the rigid rod result reported in the literature: Monolayers of rigid rods display 
neither stable tilted phases nor a liquid expanded phase. 
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We shall comment on this diagram in a few points. 

First, monolayers of chains with fixed anisotropy v , v > 16.4, display tilted 
phases only if the chains are sufficiently flexible. This substantiates our claim, 
that tilt order is stabilized by chain flexibility. 

Second, the untilted condensed phase and the expanded phase, both fluid, 
are not fundamentally different from each other. The possibility of, e.g., hex- 
atic order is ignored within our approximations. Such ordering has however 
been reported in the liquid condensed phase of lipid monolayers ||, and is 
presumably present in our model too. If this is indeed the case, the coexis- 
tence between the LE phase and the LS or a L 2 phase is expected to end in a 
multicritical point, and the transition to turn into a continuous transition at 
lower values of u or v. 

Third, the role of the temperature has to be discussed. Assuming that the 
segment density in the monolayer does not change much in the interesting 
temperature regime, the temperature enters mainly via the chain stiffness u 
and the chain anisotropy v. These parameters contain the Boltzmann factor 
1/ksT, and may have a complicated temperature dependence in addition. Let 
us neglect the latter and take u, v oc 1/T for simplicity. Under this assumption, 
(v ■ u)~ l l 2 is proportional to the temperature, and [vju) 1 ! 2 is temperature 
independent. The second quantity is interesting in it's own right, since it can 
be related to the chain length n: In a continuum approximation, where chains 
are treated as space curves of length L, L oc n and stiffness i], eta oc u with 
orientational dependent interactions V, V oc v, it can be shown that only two 
of these parameters are independent, e.g., (u/n) and (vn) f34|. Hence varying 
the chain length n has the same effect as varying (v/u) 1 ^ 2 . One can speculate 
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that this remains qualitatively true for discrete chains. 

The different possibilities for temperature dependent phase behavior can 
be read off from figure 15, which redraws the diagram of figure 14 in the axis 
variables (v/u) 1 ^ 2 and (v-u) 1 ^ 2 . For example, the phase diagrams at (v/u) 1 ^ 2 = 
3 and (v/u) 1 ^ 2 = 4 resemble figs 5 and 12, respectively. In the neighborhood 
of the fluid-fluid coexistence region, increasing the "chain length" variable 
(v/u) l l 2 produces almost the same effect than decreasing the temperature. 
This fits to the experimental observation that the addition of two (CH 2 ) groups 
to a system has a comparable effect to the reduction of the temperature by 
10 - 20°C @|. At (v/u) 1/2 = 3.4, the phase behavior of figure 13 is 



recovered, which is similar to the experimental phase diagram sketched in 
figure 1. 

Note that mean field theories generally overestimate transition tempera- 
tures. The effect is particularly strong in two dimensional systems, where the 
fluctuations even prevent the possibility of true long range tilt order (p7| ). 
and second order tilting transitions are replaced by Kosterlitz Thouless type 
transitions. Hence the phase diagrams cannot be expected to be quantitatively 
correct, and Figure 15 gives just a qualitative picture of the phase behavior. 
This picture could be tested in simulations, by systematic variations of chain 
length and chain stiffness. 

4 Conclusions. 

We have discussed the interplay of chain anisotropy and conformational en- 
tropy in simple model systems for Langmuir monolayers: Systems of short 



20 



chains, which are confined to a planar surface at one end. The phase behav- 
ior as a function of chain stiffness and effective anisotropic interaction was 
calculated in mean field approximation. 

We found that systems of chains with fixed grafting points, i.e., fixed homo- 
geneous grafting density, display tilt order in a density interval. It is bounded 
by a continuous transition to an untilted phase at high coverage, and by a dis- 
continuous transition at low coverage. The high coverage transition involves 
ordering of whole chains and is of XY type, the low coverage transition is 
caused by ordering of segments and is reminiscent of the nematic/isotropic 
transition in liquid crystals. Note that beyond mean field theory, long wave- 



length fluctuations of the direction of tilt destroy the long range tilt order f37 
However, one can still expect quasi long range order, i.e., correlation functions 
decay algebraically. 

If the chains are free to move in the plane, tilt order is replaced by phase 
separation in systems of stiff chains. In systems of flexible chains, tilted phases 
remain stable to some extent. The conformational entropy of the chains sta- 
bilizes tilt order. In fact, it favors phases at lower surface coverage in general, 
which engenders both tilted phases and an additional untilted phase, the liquid 
expanded phase. 

As a function of the chain stiffness (or, as we have argued, the chain length), 
one can distinguish between four different regimes. 

(a) Very stiff chains (rigid rod limit): Only one first order order transition is 
found, from the highly diluted gas phase to the untilted liquid condensed 
phase. 

(b) Stiff chains (or short chains): An additional untilted phase appears in a 
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temperature interval. One finds two successive fluid fluid transitions from 
the gas phase, passing the liquid expanded phase, to the liquid condensed 
phase. 

(c) Chains of intermediate stiffness: Tilted phases can be stable. Depending 
on the temperature, the liquid expanded phase coexists with either a 
tilted or an untilted condensed phase. 

(d) Flexible chains (or long chains): The liquid expanded phase coexists with 
a tilted condensed phase. Upon compression of the monolayer, the tilted 
phase turns into an untilted phase via a continuous or first order transi- 
tion. 

Hence a complex phenomenology is found already in this simple model, 
which incorporates only a few aspects of the hydrophobic tails in amphiphilic 
molecules, and entirely disregards the structure of the head groups. The dif- 
ferent phases in Langmuir monolayers at low surface coverage are largely re- 
covered. 

We conclude that the essential features of the phase behavior of Lang- 
muir monolayers can already be produced by the alkane tails of the surfactant 
molecules alone. Nevertheless, the head groups have an important influence on 
the phase diagram. For example, it has been mentioned, that tilted phases can 
be stabilized by a mismatch between head group and tail segment size. This 
is most likely the dominant tilting mechanism in monolayers of single chain 
amphiphiles, e.g., fatty acids. Future investigations will have to explore this 
possibility. 
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Figure Captions 



Figure 1: Phase diagram of Langmuir monolayers at low surface coverage 
(schematic). The liquid-gas coexistence region is represented in a com- 
pressed way relative to the liquid expanded-liquid condensed coexistence 
region. Whether the latter ends in an upper critical, or turns into a second 
order transition (indicated by the dashed line) in a multicritical point, has 
yet to be established (after Ref. [[3]). 

Figure 2: Schematic picture of the model. Inset shows functional form of the 
bending potential. 

Figure 3: Free energy per particle vs molecular area at chain stiffness u = 2 
and anisotropy v = 13.7. Thin line indicates the Maxwell construction. 
Inset shows the Gibbs free energy per area g/(/ceT) vs. the chemical 
potential. 

Figure 4: Free energy per particle vs molecular area (a) for u = 2 and different 
values of v; (b) for v = 13.7 and different values of u. In (a) different offset 
values have been subtracted from the free energy. A state with tilt order 
emerges at high chain stiffness or high anisotropy (dashed line). 

Figure 5: Phase diagram in the plane of anisotropy v and molecular area A 
at chain stiffness u = 2. 

Figure 6: Free energy per particle vs molecular area at chain stiffness u = 1 
and anisotropy v = 20.3. Two solutions of the mean field equations are 
shown, one corresponding to an untilted state (thick solid line) and one 
describing a tilted state (dashed line). Thin line indicates the Maxwell 
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construction. 

Figure 7: In-plane alignment of segments d\\ vs molecular area at w = 1, 
v = 20.3. Solid line corresponds to the untilted state, dashed line to the 
tilted state. At fixed grafting density, the tilted state is stable in the 
coverage region between I and II. 

Figure 8: Nematic order parameter S and biaxiality rj vs molecular area at 
u — 1, v — 20.3. Solid line shows results for the untilted state, dashed 
line for the tilted state. Also indicated are the locations of the tilting 
transitions I and II at fixed grafting density, and of the coexistence regions 
between liquid phases in systems of mobile chains. 

Figure 9: Types of tilting transitions (see text for explanation). 

Figure 10: Density profiles of the monolayer in the direction z perpendicular 
to the interface, at u — 1, v — 20.3 in different phases (different molecular 
areas A). Long and short dashed lines show the center of mass densities 
of tail and head segments ph,t( z ), respectively; solid line shows the total 
segment density p(z). 

Figure 11: Difference between the distribution of bending angles P(cos9) in 
the middle (i.e., between second and third tail segment) and at the end 
of the chains, AP mid (cos9) = P mi d{cos9) — P en d(cos9). Results are shown 
for the parameters u — 1, v — 20.3 and different states (stable or unstable) 
at different molecular areas A. Inset shows the distribution P en d of the 
outermost angle, which was identical in all cases. 

Figure 12: Phase diagram in the plane of anisotropy v and molecular area A 
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at chain stiffness u — 1. Inset shows a blow-up of the region where the 
liquid expanded phase is stable. 

Figure 13: Phase diagram in the plane of anisotropy v and molecular area A 
at chain stiffness u — 1.5. 

Figure 14: Projection of the phase diagram in anisotropy v, chain stiffness u 
and molecular area A into the (u, v) plane. Short dashed lines indicate 
multicritical lines, long dashed line critical lines, and the solid lines are 
triple lines, where three phases can coexist as indicated. Shaded areas 
are parameter regions where tilted phases can be stable. Coexistence of 
liquid phases is found in the hatched area (expanded phase and one of 
the condensed phases) and in the dark shaded area (tilted and untilted 
condensed phase). See text for further explanation. 

Figure 15: Same as figure 14, with different axis variables. See text for ex- 
planation. 
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